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Abstract 

We compute various kinematical distributions for one-jet and two-jet inclusive photo- 
production at HERA. Our results are accurate to next-to-leading order in QCD. We use the 
subtraction method for the cancellation of infrared singularities. We perform a thorough 
study of the reliability of QCD predictions; in particular, we consider the scale dependence 
of our results and discuss the cases when the perturbative expansion might break down. 
We also deal with the problem of the experimental definition of the pointlike and hadronic 
components of the incident photon, and briefly discuss the sensitivity of QCD predictions 
upon the input parameters of the calculation, like a s and the parton densities. 
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1. Introduction 

Jet production is a frequent phenomenon in high-energy collisions. It is character- 
ized by large rates, which allow measurements of its kinematical features even in the 
case of high jet multiplicity (like six-jet events at Tevatron and five-jet events at LEP). 
Complete next-to-leading order QCD calculations are available for one- and two-jet 
inclusive quantities in hadronic collisions |ffj, and up to four-jet inclusive quantities in 
e + e~ collisions ||. The typical energy scale for the hard process is of the order of the 
average jet transverse momentum; at presently operating colliders, this means that 
the hard production process takes place in an energy domain where the predictions 
of perturbative QCD are particularly accurate. 

At hadron colliders, for one- and two-jet inclusive observables the uncertainty 
affecting the theoretical results is smaller than experimental errors, and stringent 
tests for QCD can be carried out. Although most of the comparisons between theory 
and data are quite successful, some issues need to be better clarified. In particular, 
we recall that the CDF collaboration reports an excess of events in the tail of the 
E T distribution with respect to the theoretical predictions (this discrepancy, however, 
can be absorbed in a suitable modification of the parton densities ffl), and that the 
comparison between data || taken at different center-of-mass energies in pp collisions 
is still not consistent with the scaling found in QCD. 

The amount of experimental information relevant for lepton-hadron collisions is 
presently not as detailed as in the case of hadronic or e + e~ processes. Nevertheless, in 
the near future the increased luminosity of the HERA collider will allow a statistically 
significant study also for this kind of hard scattering. In jet production at HERA, one 
can consider two different kinematical regimes. In one class of events, the electron is 
scattered at large angles (deep inelastic scattering). In another class of events, the 
electron is lost in the beam-pipe (or tagged at very small angles); in this case, it acts 
as a source of real photons, which eventually interact with the partons in the proton 
beam. It is well known that an on-shell photon has a finite probability of fluctuating 
into a hadronic state before undergoing a hard collision. Therefore, for a comparison 
with experimental data on jet photoproduction, one must consider QCD processes in 
which the photon directly enters the hard parton scattering, and processes in which 
the photon emits partons (in a way parametrized by some universal but non-calculable 
parton densities) which enter the hard parton scattering. 
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Partial next-to-leading order QCD results for jet photoproduction have been avail- 
able for some time ||. Complete calculations for one-jet and two-jet inclusive quan- 
tities which include the contributions of all the relevant partonic subprocesses have 
been reported in refs. Comparison with data are encouraging (see for exam- 

ple refs. P!0|JII| for a recent review on this point), but larger statistics is needed in 
order to study more exclusive quantities and to fully test the predictions of the the- 
ory. Therefore, the aim of this paper is not to perform a comparison with existing 
data, but rather to discuss some general features of the theoretical predictions of jet 
photoproduction in the HERA energy range. Our predictions are obtained using the 
computer codes | recently presented in ref . |§ , and are based on the general formalism 
of ref. |TI|. The fundamental difference between the results of refs. |T],[8] and those 
presented here is that the former are obtained in the framework of the slicing method, 
while the latter rely upon the subtraction method. Both methods have been devised 
in order to cancel analytically the infrared divergencies which arise in the intermediate 
steps of any next-to-leading order QCD calculation. The slicing method requires the 
matrix elements for the real emission processes to be approximated in those regions 
of the phase space which are close to an infrared (soft or collinear) singularity. These 
regions are defined by means of some non-physical parameters. Physical results must 
be independent of these parameters; it is in general a non-trivial task to prove (nu- 
merically) that this is indeed the case, at the required level of accuracy. Notice also 
that in principle this proof has to be carried out for each observable, since different 
quantities may display different convergence properties. In the subtraction method, 
the exact expression of the matrix elements is used in the whole phase space. The 
main advantage of this method is that it does not require the introduction of any non- 
physical parameter. In the following, we will encounter cases where this feature of the 
subtraction method will turn out to be advantageous in the calculation of physical 
quantities. 

This paper is organized as follows. In section 2 we study the problem of the re- 
liability of the next-to-leading order QCD predictions for jet observables at HERA, 
where the typical transverse momenta are sizeably smaller than in the case of hadron 
collider experiments. We consider several one-jet and two-jet inclusive quantities, 
and we isolate the regions where QCD perturbative expansion breaks down, and an 
all-order resummation of soft-gluon effects is mandatory. We also discuss the prob- 

2 The codes, which return parton kinematical configurations with an appropriate weight, are 
available upon request. 
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lem of the operational separation of the pointlike (or direct) and of the hadronic (or 
resolved) components of the cross section. In section 3 we present some predictions 
for quantities whose measurement would help the understanding of the jet produc- 
tion mechanism and ultimately to better test QCD in the context of photon-hadron 
collisions. Finally, in section 4 we report our conclusions. 



2. General features of QCD predictions 

Any production cross section in electron-proton collisions is dominated by the 
exchange of low-virtuality photons. The electron behaves therefore as a broad-band 
beam of real photons, whose momenta are distributed according to the Weizsacker- 



Williams function [13] 
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where m e is the electron mass and y is the fraction of the electron longitudinal mo- 
mentum carried by the photon. Notice that in eq. ( |2.1| ) terms which are non-singular 
for m e — > cannot be neglected in the HERA energy range, because of the 1/y singu- 



larity [|i4| , |15| . The quantity Q must be taken to be the minimum between the typical 
energy scale which characterizes the production process, and the upper limit of the 
absolute value of the photon virtuality (see ref. $M for a detailed discussion of this 



point). The cross section for a generic electron-proton scattering process is given by 
da ep (K e , K p ) = f VmaX dy f<f\y) da lp {yK e , K p ), (2.2) 

where K e and K p are the momenta of the incoming electron and proton respectively, 
and < y min < y-max < 1 are fixed by kinematical boundary conditions or experi- 
mental cuts. The relevant hard scattering quantity is therefore the cross section 
for photon-proton collisions, da lv . Since the HERA experiments can measure the 
momentum of the scattered electron down to very small values of the scattering 
angle, it is possible, in principle, to measure directly da lp for a fixed value of the 
photon energy (with the current electron and proton energies, E e = 27.5 GeV and 
E p = 820 GeV, the accessible range in the photon-proton center-of-mass energy is at 
best 100 GeV< E cm (^p) < 280 GeV). 

The possibility of measuring the photon-proton jet cross section for several differ- 
ent center-of-mass energies would be an interesting test of the scaling laws predicted 



-4- 




Figure 1: Comparison between monochromatic photon-proton and electron- 
proton (Weizsacker- Williams approximation) cross sections, for various jet 
distributions. 

by QCD. For this reason, we have studied the energy dependence of differential jet 
distributions in monochromatic photon-proton collisions, comparing them with the 
same quantities obtained in electron-proton collisions in the Weizsacker- Williams ap- 
proximation, eq. (|2.2| ). The results, obtained at NLO in QCD, are shown in fig. [1]. In 
the two upper plots we present single-inclusive distributions in the transverse energy 
E T of the jet and in its pseudorapidity r\ for E T > E^ ut , E^ ut = 10 GeV. In the two 
lower plots, we show the distributions in the transverse momentum, pij , and invariant 
mass, Mjj, of the pair of the two hardest jets in the event. In this case, we require 
the transverse energies of the two observed jets to satisfy the constraints E\ T > E™*, 
E 2t > E^, with E^ = E$f = 10 GeV (the fact that the same value of the transverse 
energy cut is used for both jets will be discussed below in detail). Jets are defined by 
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the cone algorithm 0] with opening angle R — 1. We used the sets MRS A' [17| and 
GRV-HO ]18[ for parton densities in the proton and the photon respectively. The 
renormalization and factorization scales have been chosen equal to half the trans- 
verse energy of the event; from now on, this choice will be our default and will be 
denoted as hq. A qcd has been set equal to the value suggested by the MRS A' parton 
densities, A^ s = 152 MeV. Each distribution in fig. [I] is presented for two different 
energies of the photon-proton system, 134 GeV (dotted curves) and 268 GeV (dashed 
curves), which correspond to y = 0.2 and y = 0.8 in eq. (|2.2|) respectively. The solid 
curves show the predictions for electron-proton collisions at E cm (ep) = 300 GeV in 
the Weizsacker- Williams approximation, eq. ( |2.2|) , with y min = 0.2, y ma x = 0.8, and 
Q 2 = 4 GeV 2 , rescaled by a factor 
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26.9. (2.3) 



The difference in normalization between the curves relevant for the two photon- 
proton center-of-mass energies is sizeable. Phase-space effects are clearly visible in 
the tail of the E T , p>J and Mjj distributions, where the prediction for E cm = 134 GeV 
decreases faster than the power-like behaviour typical of QCD, and in the range of 
negative i] (the photon is coming from the right, as in the HERA conventions). By 
construction, there is no dynamical information on the hard scattering process in 
the Weizsacker- Williams approximation, which is simply a weighted sum of 7p cross 
sections; this is also clear from the shape of the solid curves in fig. |l| For this reason, 
in this section we will mainly consider monochromatic photon-proton collisions at 
E cm = 268 GeV. We have verified that our conclusions apply down to the lowest 
center-of-mass energies relevant at HERA. 

We now turn to the problem of studying the reliability of the perturbative expan- 
sion. To this purpose, we have considered photon-proton collisions at 268 GeV as a 
representative case and we have computed various distributions, varying factorization 
and renormalization scales by a factor of 2 around our default choice, as customary in 
QCD in order to assess the magnitude of the neglected higher-order corrections. QCD 
predictions for photon-hadron cross sections are computed in terms of two distinct 
components: the so-called pointlike component, in which the photon directly inter- 
acts with partons in the hadron, and the hadronic component, in which the photon 
fluctuates into hadronic states, which subsequently interact with the hadron. This 
distinction is well defined at leading order, but becomes arbitrary when higher order 



-6- 



terms in the perturbative expansion are included (we will thoroughly discuss this 
problem in the following). Explicitly, the cross section is written as 

da, p (K„ K p ) = daV° p int (K„ K p ) + da^ dl (K„ K p ) (2.4) 

where 

1 ftw ft HI ftl ), (2-5) 

3 

1 ftw ft HI ftj) 1 

(2.6) 

and fi (fi ) are the distribution functions of parton i inside the proton (photon). 
da 1 j and da^ are the subtracted short-distance partonic cross sections, which are 
finite at any order in perturbative QCD. The dependence of da^° mt on (jl' r ,(jl' h , and 
of da^p dr on /i^, fi' H , cancels in each component separately, up to NNLO terms. On 
the other hand, the dependence on /i 7 in da^ p mt is compensated by the analogous 
dependence in da^ T (see for example ref. |0| for a detailed discussion on this point). 



In the language of Altarelli-Parisi equations, the collinear singularities of the photon 
leg in the pointlike component are re-absorbed (at a scale /x 7 ) in the parton densities 
of the photon entering the hadronic component. 

We begin by studying single-inclusive quantities and two-jet observables, defined 
with E^f 7^ E^f. We set all the scales in eqs. Q2.5|) and (|2.6f ) to the same value, which 
we denote by /i, and we study the dependence of the physical quantities upon In 
this case, this procedure gives a good estimate of the full uncertainty affecting the 
results, which in principle should be obtained by considering independent variations of 
all the scales. In the left plot of fig. ^we present the distribution in the pseudorapidity 
of jets with E T > E^ ut , for E^ ut = 10 GeV and 20 GeV. In the right plot of the same 
figure we show the distribution in the invariant mass of the two hardest jets, for E™* = 
10 GeV, E™ 1 = 15 GeV (upper curves) and E%* = 20 GeV, E^ = 25 GeV (lower 
curves). The jets are defined by the cone algorithm with R = 1. The scale dependence 
is moderate, and slightly reduces with increasing transverse energy cuts. We checked 
that other kinematical quantities, like the E T of the jet, the pseudorapidity of the 
pair, the difference in rapidity of the pair, show a similar scale dependence. From 
the left-hand side of fig. we see that the scale dependence of the pseudorapidity 
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Figure 2: ^ca/e dependence of the single-inclusive jet pseudorapidity and two- 
jet invariant mass distributions, for different values of the minimum allowed 
transverse energy of the observed jets. 
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Figure 3: Scale dependence of the azimuthal distance distribution in two-jet 
events, for different values of the minimum allowed transverse energy of the 
two hardest jets (E^f / Egj?). 
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distribution is slightly stronger for 77 values which are far from the central region. 
This is what we expect, since large 77 values correspond to small transverse energies, 
and therefore to less reliable QCD predictions. 

We also computed the same quantities for a cone with an opening angle of R = 0.7. 
In this case, the scale dependence of both the pseudorapidity and the invariant mass 
distributions is slightly reduced. 

The situation is different for the azimuthal correlation A<f>jj between the two jets 
with largest transverse momenta, shown in fig. |J| For Afyj < n the two sets of 
curves corresponding to the different choices of transverse momentum cuts show the 
same scale dependence, larger than in the case of pseudorapidity and invariant mass 
distributions of fig. [2| This is due to the fact that the A<ftjj correlation is a pure NLO 
effect in this region. The scale dependence reduces for Afyj ~ n, as can be also seen 
from the small inserted figure. 

The A(j)jj correlation computed with a cone of R = 0.7 shows the same scale 
dependence as in the case R = 1 for A<pjj < 71", and a stronger scale dependence for 

Afaj ~ 7T. 

We now consider the problem of computing two-jet inclusive quantities for Ef!f = 
= E™ f . As far as infrared safeness is concerned, there is nothing special in 
this choice. The cross section is well-defined and finite at any order in perturbation 
theory. On the other hand, there are quantities which at next-to-leading order display 
a pathological behaviour. This can be seen very easily by studying the inclusive two- 
jet total cross section 

(7 2 (A) = a(E lT > Ef, E 2T > E* ut + A) (2.7) 

for A — > 0. The next-to- leading order QCD results are shown in fig. |], where we have 
chosen Eg* = 10 GeV and Eg* = 20 GeV (the curve for Eg* = 20 GeV has been 
rescaled by a factor of 11 in order to make both curves visible on the same plot). 
Notice that the value of 02 (A) for A = is finite, as expected for an infrared-safe 
quantity. Observe also, on the other hand, that (72(A) has an infinite slope in A = 0. 
This fact can be understood in the following way. We consider the real emission 
contribution when one of the emitted partons is quasi-collinear to one of the initial 
state partons (at NLO in QCD this implies that the two jets are identified with the 
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Figure 4: Inclusive two-jet total cross section for E\ T > E^ and 
> E™ 1 + A, as a function of A, for two different values of ' E^ 1 . 

two other partons). The leading collinear singularity of this contribution is given by 
a<r) = [ d 2 PlT 6(E 1T - [ d 2 p 2T 9(E 2T - E^ - A) \ (2.8) 

J J \PlT + P2t\ +0 

where 5 acts as a collinear cutoff. The integral can be easily computed, and we get 

a (r) = A(A,S) + B\og5-C- (A + 5) log(A + <5), (2.9) 

where both A(A, 6) and its first derivative with respect to A are regular in A = 
for any 5, including 5 = 0, and C is a positive coefficient. The term B log 5 is the 
genuine collinear singularity, and it is cancelled by the corresponding singular terms 
in the two-body contribution. For fixed A / 0, one can safely take 5 = in the 
last term. For A —>■ this term vanishes, but its first derivative with respect to 
A diverges, and this is the origin of the behaviour observed in the full calculation, 
shown in fig. Alternatively, one could have started with A = 0, as we did for the 
distributions shown in fig. |l|. In this case, the last term in eq. ([2.9D becomes 5 log 5, 
which vanishes for 5^-0, although less rapidly than terms linear in S. Because of this 
extra S dependence, the A = case is a reason of concern || when the calculation is 
performed with a technique which requires keeping 5 ^ 0, as for example the slicing 
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method. We remark that in the subtraction method the introduction of the non- 
physical parameter S is avoided, and therefore the A = case does not require any 
special care. 

One striking feature of fig. f| is that 02(A) is not a monotonically decreasing 
function for increasing A, as one might expect on the basis of simple phase-space 



considerations. This behaviour is also present is our simple model, eq. (2J)), because 
of the negative sign in front of last term. This effect is induced by the truncation of 
the perturbative series at NLO. Roughly speaking, at this order there is not enough 
emission of soft real gluons to compensate for the large negative contribution of the 
soft- virtual terms. We therefore expect quantities like 02(A) to be poorly predicted by 
NLO QCD for A = 0. However, the mismatch between virtual and real contributions 
is effective only in some special regions of the phase space; typical examples are the 
threshold of the invariant mass of the pair and the region around A<pjj = it in the 
azimuthal distance correlation of the two hardest jets. It is well known that fixed- 
order perturbation theory is not reliable in these regions, and that resummation of 
all-order contribution is needed. Despite this fact, the cross section is expected to be 
well-behaved in the remaining part of the phase space. 

In the following, we will show that two-jet inclusive distributions can be safely 
studied even in the case of equal transverse energy cuts on both jets. To address this 
issue, we performed a more careful study of the scale dependence of the results. In 
particular, we varied independently /x 7 and fi R = fi H , where fi R = fjf R and \i H = fi' H 
when we consider the pointlike component, eq. Q2.5Q , while fi R = jj!' R and fi H = n" H 
when we consider the hadronic component, eq. ( |2.6| ). We also verified that indepen- 
dent variations of fi R and \i H would not modify our conclusions. We begin with the 
A(pjj correlation, shown in fig. |5|. The region A<pjj < n does not pose any problem, 
being a pure NLO effect. For this reason, the curves were calculated again with 
jjL^i —— j^ R l^H' As in the case E™* 7^ E!^, no improvement in the scale dependence 
is seen when transverse energy cuts are increased. The small inserted figure shows the 
region A<f>jj ~ 7r; in this case, the curves have been evaluated by varying the scales in- 
dependently. In the bin around A<pjj = ir the curves corresponding to E^ ut = 10 GeV 
show a very strong scale dependence, thus signalling a failure of the perturbative 
expansion in this region. Notice also that most of the uncertainty band in the last 
bin gives negative cross sections. By considering independent scale variations, this 
band increases of about 25% with respect to the case when all the scales are set to 
the same value. When E^ ut is increased to 20 GeV, this effect is more moderate, but 
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Figure 5: Scale dependence of the azimuthal distance distribution in two- jet 
events, for different values of the minimum allowed transverse energy of the 
two hardest jets (E{f = E$? = E° ut ). 

the scale dependence is still larger than in the region A<f>jj < it (this is consistent 
with fig. |], where we see that the total cross section decreases faster for A — > in 
the case of smaller transverse momentum cuts, thus indicating that in this case soft 
gluon emission is more important). We expect that an all-order resummation would 
have a moderate impact in this case. Of course, the above considerations are partially 
dependent on the chosen bin size: enlarging the bin around A<pjj = it enough, the 
perturbative prediction becomes reliable for any choice of E^ ut . 

We now turn to the distribution of the invariant mass of the pair. This distribution 
is non-trivial already at leading order, and therefore, in principle, the effect of soft 
gluon emission is not confined in some particular region of the phase space. To clearly 
show the effect of an independent variation of the scales, we separately present in fig. |6] 
the pointlike and the hadronic contributions, for the choice = E^f = 10 GeV. 
The solid curves are the defaults, the dotted (dashed) curves are obtained by varying 
H>r = A*h (/^y) • We see that the pointlike and the hadronic component display opposite 
behaviours with respect to the scale /i 7 (for /i 7 = 2/i , the curve is larger than the 
default one for the hadronic component, and smaller for the pointlike component). 
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Figure 6: Scale dependence of the invariant mass distribution in two- jet 
events. The pointlike (left) and hadronic (right) components are separately 
shown (E{f = E^ f = E^ ut ). 



This shows the cancellation mechanism between eqs. ( |2.5|) and fl2.6|) . In particular, 
the fact that the hadronic component is larger than the pointlike one implies that a 
moderate /i 7 dependence in the hadronic contribution induces a huge /i 7 dependence 
in the pointlike contribution, as can be seen from fig. []. This does not mean that the 
prediction of the theory is not reliable: in fact, the only measurable quantity is the sum 
of the two components, shown in fig. |7|. From that figure, we see that the compensation 
between the /i 7 dependence of the two components is remarkable, and the residual 
dependence is rather small. On the other hand, there is no cancellation mechanism 
in the case of the fi a and \i H dependence, since these quantities are different in the 
pointlike and hadronic components (see eqs. ( |2.5|) and Q2.6Q ). The corresponding 
uncertainties must be summed incoherently. Nevertheless, these uncertainties are 
small, therefore resulting in a small uncertainty on the physical result, fig. |7|. 

This discussion does not apply to the region where Mjj is close to the threshold. 
From fig. [7|we see that the bin around Mjj = 20 GeV shows a large scale dependence, 
and that some of the curves are negative. In this region, the perturbative expansion 
therefore breaks down, like in the case A<ftjj ~ ir we discussed above. 

We conclude that NLO QCD provides reasonably stable predictions for jet produc- 
tion at HERA. Single-inclusive variables are very well described. In the case of two-jet 
inclusive quantities, we have to distinguish the cases ^ E^ and E^f = E^. In 
the former case, the perturbative result is rather accurate, and gives reliable predic- 
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Figure 7: Scale dependence of the invariant mass distribution in two-jet 
events (E$* = E%* = E^ ut ). 



tions in the full range of any quantity (this assumes that E^"* and are not too 
close to each other. By inspection of fig. f|, it is safe to take \E%£ — E%£\ > 2 GeV). In 
the latter case, there are regions of the phase space which are particularly sensitive to 
soft gluon emission, and therefore require an all-order resummation. In these regions 
(A0jj ~ 7r, p>J ~ 0, Mjj close to the threshold) the fixed-order QCD results are not 
reliable, and no significant comparison can be made with experimental results (which, 
on the other hand, can be safely obtained, since E^ = E$f does not imply infrared 
non-safeness) . Elsewhere, two-jet correlations can be predicted to a level of accuracy 
comparable to the one obtained in the case E^f ^ E^ 

As discussed above (see in particular fig. the pointlike and hadronic components 
are strongly correlated in perturbative QCD beyond leading order. This issue is inter- 
esting in the light of the fact that some analyses of experimental data are performed 
with specific kinematical cuts, imposed in order to measure physical quantities which 
are supposed to be dominated by the pointlike or the hadronic component of the cross 
section for any reasonable scale choice. In jet physics, such cuts are usually defined 
in terms of the variable 

= 1T _ 2T , (2.10) 



2E~, 



-14- 




Figure 8: Scale dependence of the invariant mass distribution, for x 7 > 
0.75 (left) and x 7 < 0.75 (right), compared with the pointlike and hadronic 
component respectively. 



where E iT and rji (i = 1,2) are the transverse energies and pseudorapidities of the 
two hardest jets in the event. For two-jet production at leading order, x 7 is equal to 
1 for the pointlike component (in the hadronic component x 1 is the fraction of the 
photon longitudinal momentum carried by the interacting parton). At higher orders, 
this identification no longer holds, but one might expect that the high-x 7 region 
is dominated by the pointlike component, and the low-x 7 region by the hadronic 
one. Following this criterion, in some analyses the pointlike (hadronic) component is 
operationally defined as the measured cross section for x 7 larger (smaller) than some 
fixed value, usually 0.75. With such a definition, one might for example hope to pose 
direct constraints on the parton densities in the photon, or to get information on the 
details of the underlying parton dynamics, using jet data. 

In order to assess the effectiveness of this procedure, we have computed the pair 
invariant mass distribution imposing the cuts x 1 > 0.75 or x 7 < 0.75. We considered 
again 7J9 collisions at a center-of-mass energy of 268 GeV, defining the jet by the cone 
algorithm with R = 1. In fig. § we present our results for the choice = 10 GeV, 
= 15 GeV. In the figure on the left (right), we show our prediction for the curve 
with x 1 > 0.75 (s 7 < 0.75), together with the pointlike (hadronic) component. As 
before in the case E™* ^ E^, the scale dependence has been studied by setting all 
the scales to the same value. The cuts — 1 < 771,772 < 2 have been imposed on the 
pseudorapidity of the two observed jets in order to simulate a realistic geometrical 
acceptance. We observe that the hadronic component (see fig. |8] right) is sizeably 
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Figure 9: Pseudorapidity distribution for x 7 > 0.75 (left) and x 7 < 0.75 
(right), compared with the pointlike and hadronic component respectively. The 
scales have been set to the default value. 



larger than the distribution for x 7 < 0.75. This implies that in the region x 7 < 0.75 
the pointlike and the hadronic contributions mix significantly for any scale choice. 
The same pattern is displayed by the pointlike component and the distribution for 
x 7 > 0.75 (fig. ||left); in this case the scale dependence of the pointlike component 
is somewhat larger. We also computed the distributions of fig. || in the case = 
E^p = 10 GeV. We observe the same behaviour; the difference between the hadronic 
(pointlike) component and the curve for x 7 < 0.75 (x 7 > 0.75) is even larger than 
that shown in fig. [8]. 

We found that the mixing between the hadronic and pointlike components is 
sizeable for all the physical observables, even in simple cases like single-inclusive dis- 
tributions. As an example, we present in fig. |9] the comparison between the hadronic 
(pointlike) component of the rj distribution, and the corresponding curve obtained 
with ai 7 < 0.75 (x 1 > 0.75) cut. Although in this case the shapes are quite similar, 
the absolute normalization is rather different. We have checked that the difference 
between the two curves in each plot of fig. |9| is larger than the corresponding scale 
uncertainty. 

This discussion suggests that, because of the effect of radiative corrections, a cut 
at x 7 = 0.75 is probably not very useful to distinguish between partonic subprocesses 
with or without an incoming photon. Nevertheless, it is still conceivable that by 
cutting in x~, one may select observables which are particularly sensitive to the parton 
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Figure 10: Dependence of the invariant mass distribution upon the parton 
densities in the photon. 



densities in the photon. In order to further investigate this issue, we have computed 
the invariant mass distribution for x 7 < 0.75 with rather an extreme choice of the 



parametrization of parton densities in the photon, namely the LAC1 |20[ set, and 
compared it to our previous results obtained with the GRV-HO set. The invariant 
mass is particularly suitable to this purpose, since it is directly related to the Bjorken- 
x region probed in the collision (for a given center-of-mass energy, small invariant 
masses correspond to the small- x region). One can see from fig. [10] that the difference 
between the results obtained with the GRV and LAC1 sets is comparable to the 
uncertainty band induced by the scale variations we discussed above. In the central 
pseudorapidity region the sensitivity to the different parametrizations is very small. 
In particular, the right plot of fig. [II] shows that in this rj region the x 7 cut does 
not enhance the dependence upon the photon parton densities. Notice that fig. [1^ is 
relevant for monochromatic photon-proton collisions at E cm = 268 GeV; for smaller 
center-of-mass energies or for electron-proton collisions in the Weizsacker- Williams 
approximation the sensitivity is even more suppressed, because larger values of the 
Bjorken-x are probed. The main difference between the GRV and LAC1 sets is in the 
behaviour of the gluon density at small x. To effectively probe this region at HERA 
smaller values of the invariant mass must be chosen, and QCD predictions become 
less reliable. 
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yp collisions, E cm =200 GeV, E T °"'=6 GeV 

dashed: x,<0.75, M„>23 GeV, |?)|<0.5 
solid: x 7 >0.75, M,j>23 GeV, |»)|<0.6 
+ : x ? <0,?5, no additional cuts 
□: x^>0.75, no additional cuts 
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yp collisions, E cnl = 26B GeV, £,""'=6 GeV 

dashed: x 7 <0.75, M„>23 GeV, |j)|<0.6 
solid: x 7 >0.75, Mj,>23 GeV, [tj |<0.5 
+ : x ? <0,75, no additional cuts 
□: x r >0.75, no additional cuts 




Figure 11: Distribution in cos6*, for ^yp collisions and various kinematical 
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3. Phenomenology of jet production at HERA 

In this Section we will discuss the predictions of perturbative QCD for some 
quantities of phenomenological interest. We begin by considering the distribution in 
the variable cos 8*, defined as 



cos 6* = tanh 



(3-1) 



where r]i j2 are, as before, the pseudorapidities of the two hardest jets. In QCD at 
leading order the hadronic component is predicted to behave as (1 — |cos#*|) -2 when 
cos 6* —>■ 1, as opposed to the (1 — |cos ^* I) -1 behaviour of the pointlike component. It 
is interesting to check whether the measured distributions for x 1 < 0.75 and x 7 > 0.75 



display a similar behaviour. In fig. [11] we present the cos 9* distribution in collisions 
for E™ 1 = = 6 GeV, for two different center-of-mass energies. We see that for 
both values of E cm the curves for x 7 > 0.75 and x 1 < 0.75, represented in fig. [LI] by 
squares and crosses respectively, have approximately the same shape (the two curves 
have been normalized to have the same integral in the first four bins, in order to 
compare the shapes). This confirms the expectations of the previous section, where 
we have shown that distributions in the x 7 < 0.75 (or a; 7 > 0.75) region are the 
result of a significant mixing between the hadronic and the pointlike components. 
It is interesting to notice that a difference in the behaviour at cos#* ~ 1 of the 
two distributions obtained for x 7 < 0.75 and x 1 > 0.75 can be induced by other 
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kinematical cuts. To show this, we have computed the cos 8* distribution in the two 
x 7 regions, applying the cuts Mjj > 23 GeV and \ fj\ < 0.5, with fj = (771+772) /2. These 
are the same cuts adopted by the ZEUS collaboration in the analysis of ref. |[21|| . This 
curves (normalized as before to the same integral in the first four bins) are also shown 
in fig. |TT|; we see that the curves for x 7 < 0.75 (dashed) and x 7 > 0.75 (solid) display 
a different behaviour for cos#* — > 1; this difference is generated by the introduction 
of the kinematical cuts. This effect can be understood observing that 

Mf d = 2E lT E 2T [cosh(r/ 1 - r/ 2 ) - cos Afaj] . (3.2) 

Large values of the invariant mass correspond therefore to large values of cosh (7/1 — 7/ 2 ); 
on the other hand, 

1 + cos 2 9* . . 

cosh (771 - 772) = 3.3 

1 — cos^ V* 

It is clear from eq. ( |3.3| ) that a large- Mjj cut enhances the region cos 8* ~ 1. This 



shows that the difference in shape between the solid and dashed curves in fig. [TT 
is not directly related to a different dynamical production mechanism. The curve 
corresponding to x 7 < 0.75 has a steeper rise basically because the invariant mass 
distribution relevant for this x 7 region is softer than the invariant mass distribution 
for x 7 > 0.75. One may then argue that this is still a signal of different production 
mechanisms. But, as we discussed in the previous section, the mixing between the 
pointlike and the hadronic components is especially relevant in the case of — E^^ : , 
and therefore no particularly significant information on the parton dynamics can be 
extracted from the invariant mass distribution with a x 7 cut. 

For a given observable, the largest statistics will be collected by integrating over 
the energies of the incident photons. For this reason, we will now present predictions 
obtained for jet production in ep collisions in the Weizsacker- Williams approxima- 
tion, with E cm {ep) = 300 GeV, Q 2 = 0.01 GeV 2 and 0.2 < y < 0.8 (see eqs. Q 
and ( p.2|) ). In order to perform a realistic study, we will apply kinematical cuts which 
approximately reproduce the experimental conditions of the HERA experiments. The 
pseudorapidity of the observed jets will be restricted in the range — 1 < 77 < 2. In 
the case of two-jet quantities, we will require E^ = 11 GeV and E^ = 14 GeV. By 
definition, the two jets are always the hardest of the event. 

As discussed in the previous section, with this choice of cuts on transverse energies 
we expect the perturbative expansion to be well-behaved over the whole phase space. 
For this reason, we study the scale dependence of our results by setting all the scales to 
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Figure 12: Single-inclusive transverse energy and pseudorapidity. Jets are 
defined by the cone algorithm (solid curves; the scale dependence is also shown) 
and by the prescription of ref. (circles). 



the same value \x. We present predictions obtained with two different jet definitions, 
namely the cone algorithm, with R — 1, and the algorithm proposed by Ellis and 
Soper in ref. [22|, with D — 1. In the case of the cone algorithm we show the 
full scale variation computed as before, while for the Ellis-Soper algorithm we only 
present the predictions corresponding to our central scale choice. We recall that 
a modification of the standard cone algorithm has been introduced in ref. [23] to 
improve the comparison with data; this modification however would not change our 
conclusions, and we will use R sep = 2R in the following. 



In fig. [12] we show single-inclusive distributions in transverse energy and pseudora- 
pidity of the jet (in the latter case, we require the jet to have transverse energy larger 
than 14 GeV). In fig. Oand fig. |TJ]we show two-jet correlations which are non-trivial 
at leading order: the invariant mass of the pair, the absolute value of the difference 
in pseudorapidity of the two observed jets, and the average pseudorapidity of the 
pair, fj = {rji + ^ 2 )/2; notice that the latter quantity is equal, at leading order, to the 
rapidity of the pair. Finally, fig. [15| shows correlations trivial at leading order, the 
transverse momentum of the pair and the azimuthal distance. This set of observables 
gives a fairly complete description of the production mechanism; it is clearly possible 
to consider even more exclusive quantities, by imposing additional kinematical cuts. 

The results presented in figs. [12|, [H| and [14] display a remarkable stability with 
respect to scale choice. The difference between the default curve (solid histogram) 
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Figure 13: Invariant mass of the pair of jets with the largest transverse 
momenta. Jets are defined by the cone algorithm (solid curve; the scale de- 
pendence is also shown) and by the prescription of ref. Jgl[ / (circles). 

and the curves corresponding to \i = 2/i and \i = /io/2 (dashed and dotted his- 
tograms respectively) is about 10%. We remark that the non-logarithmic term in the 
Weizsacker- Williams function, eq. ( |2.1| ), gives a negative contribution of the order 
of 7% with Q 2 = 0.01 GeV 2 , and of 5% with Q 2 = 4 GeV 2 , and is therefore non- 
negligible. Using the jet definition of ref. [p2]| , represented by circles, one gets shapes 
similar to the ones given by the cone algorithm with R = 0.7. 

Consistently with what we found in the previous section, the scale dependence is 
larger in the case of the plots of fig. which are a pure next-to-leading order effect 
except for the bins around jpj = and A<pjj = tt. In this case, differences of about 
25% are induced by scale variations. The distributions obtained with the Ellis-Soper 
algorithm are broader than those obtained with the cone algorithm and R — 1, when 
the jets tend to be more back-to-back in the transverse plane. 

We have also checked that the choice of parton densities does not influence the 
QCD predictions very much. For example, using the MRS125 set, characterized by 
the fact that a s {mz) is fixed to the value 0.125, the shape of the distributions does 
not change, and the modification in the overall normalization is easily traced back to 
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Figure 14: Average pseudorapidity and difference in pseudorapidity of the 
pair of jets with the largest transverse momenta. Jets are defined by the cone 
algorithm (solid curves; the scale dependence is also shown) and by the pre- 
scription of ref. (circles). 




Figure 15: Transverse momentum and azimuthal correlation of the pair of 
jets with the largest transverse momenta. Jets are defined by the cone algo- 
rithm (solid curves; the scale dependence is also shown) and by the prescription 
of ref. (circles). 



the different value of a s . 
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4. Conclusions 

In this paper we have studied the phenomenology of jet photoproduction at HERA 
at next-to-leading order in QCD. We used the subtraction method, which does not 
require an approximate expression for the matrix elements and the introduction of 
non-physical infrared cutoffs which are necessary in other approaches (like the slic- 
ing method). We found that, starting from transverse energy scales of the order of 
10 GeV, the theoretical predictions are sufficiently stable with respect to scale vari- 
ation to allow for a significant comparison with high-statistics data. In some cases, 
like two-jet inclusive quantities defined with equal cuts on the transverse energy of 
the observed jets, the next-to-leading order result turns out to be unreliable in some 
special regions of the phase space, corresponding typically to configurations in which 
the jets are back-to-back in the transverse plane, or to the threshold of the invariant 
mass distribution. In these cases, an all-order resummation would be needed to get 
a consistent result. On the other hand, in the remaining part of the phase space 
the next-to-leading order result is well-behaved, and meaningful comparisons with 
experimental data can be carried out. 

We have considered the problem of the separation of the pointlike and the hadronic 
components of the cross sections. This task could in principle be performed, if a suit- 
able operational definition could be given for the two components (the pointlike and 
hadronic components defined as in eqs. fl2.5|) and ( |2.6| ) are not physical observables 
beyond leading order). We found that by cutting in the variable x 7 , defined with 
the two hardest jets of a given event, the operationally defined pointlike or hadronic 
components are actually a sizeable mixture between the next-to-leading order point- 
like and hadronic cross sections, for any reasonable choice of scales. This implies that 
this kind of cut is not very useful to extract information on the underlying parton 
dynamics. Also, at the transverse energy scales where the perturbative expansions is 
reliable, the sensitivity of the operationally defined hadronic component to the parton 
densities in the photon is rather limited. 

Finally, we have considered a sample set of one-jet and two-jet inclusive observ- 
ables which can be measured in order to have a detailed description of the production 
mechanism. We have computed these quantities in the Weizsacker- Williams approxi- 
mation, applying a set of realistic cuts in order to reproduce the experimental analyses 
of the HERA experiments. We found that the results are remarkably stable with re- 
spect to scale variation, and therefore that the measurements of these quantities could 
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be used as a quality test of perturbative QCD. 
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